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^p H , Abstract. We show how non-radial motions, originating in the outskirts of clusters of 

O ■ 

!— i ' galaxies, may reduce the discrepancy between the Cold Dark Matter (CDM) predicted 

, X-ray temperature distribution function of clusters of galaxies and the observed one and 

k*" ■ also the discrepancy between the CDM predicted two-point correlation function of clus- 

' ters of galaxies and that observed. We compare Edge et al. (1990) and Henry & Arnaud 



(1991) data with the distribution function of X-ray temperature, calculated using Press- 
Schcchter's (1974 - hereafter PS) theory and Evrard's (1990) prescriptions for the mass- 
temperature relation and taking account of the non-radial motions originating from the 
gravitational interaction of the quadrupole moment of the protocluster with the tidal 
field of the matter of the neighboring protostructures. We find that the model produces 
a reasonable clusters temperature distribution. We compare the two-point cluster corre- 
lation function which takes account of the non-radial motions both with that obtained 
by Sutherland & Efstathiou (1991), from the analysis of Huchra's et al. (1990) deep 
rcdshift survey, and with the data points for the Automatic Plate Measuring (APM) 
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clusters, computed by Efstathiou et al. (1992a), showing that non-radial motions reduce 
the discrepancy between the theoretical and the obsservational correlation function. 

Key words: cosmology: theory-large scale structure of Universe - galaxies: formation 



1. Introduction 

Although at the begining the standard form of CDM was very successful in describing 
the structures observed in the Universe (galaxy clustering statistics, structure formation 
epochs, peculiar velocity flows) (Peebles, 1982; Blumenthal et al. 1984; Bardeen et al. 
1986 - hereafter BBKS; White et al. 1987; Frenk et al. 1988; Efstathiou 1990) recent 
measurements have shown several deficiencies in the model, at least, when any bias of the 
distribution of galaxies relative to the mass is constant with scale (see Babul & White 
1991; Bower 1993; Del Popolo & Gambera 1998a, 1998b). Some of the most difficult 
problems that must be reconciled with the theory are: 

the magnitude of the dipole of the angular distribution of optically selected galaxies 
(Lahav et al. 1988; Kaiser & Lahav 1989); 

the possible observations of clusters of galaxies with high velocity dispersion at z > 0.5 
(Evrard 1989); 

the strong clustering of rich clusters of galaxies, £ cc (r) ~ (r /25h~ 1 Mpc)~ 2 , far in excess 
of CDM predictions (Bahcall & Soneira 1983); 

the X-ray temperature distribution function of clusters, over-producing the observed 
clusters abundances (Bartlett & Silk 1993); 

the conflict between the normalization of the spectrum of the perturbation which is 
required by different types of observations; 

the incorrect scale dependence of the galaxy correlation function, £(r), on scales 10 to 
100 h~ 1 Mpc, having £(r) too little power on the large scales compared to the power on 
smaller scales (Maddox et al. 1990a; Saunders et al. 1991; Lahav et al. 1989; Peacock 
1991; Peacock & Nicholson 1991). 

These discrepancies between the theoretical predictions of the CDM model and the obser- 
vations led many authors to conclude that the shape of the CDM spectrum is incorrect 
and to search alternative models (Peebles 1984; Efstathiou et al. 1990a; Turner 1991; 
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White et al. 1993a, 1993b; Shafi & Stecker 1984; Valdarnini & Bonometto 1985; Holtz- 
man 1989; Schaefer 1991; Shaefcr & Shafi 1993; Holtzman & Primack 1993; Ccn ct al. 
1992; Bower et al. 1993). 

In this paper we address two of the quoted and most serious problems of the CDM model, 
namely that of the discrepancy between the predicted and observed X-ray temperature 
distribution function of clusters and that of the discrepancy between the CDM predicted 
two-point correlation function of clusters of galaxies and that observed, showing how 
them may be reduced when non-radial motions, that develop during the collapse process, 
are taken into account. 

X-ray studies of clusters of galaxies have provided a great number of quantitative data 
for the study of cosmology. The mass of a rich cluster is approximately lO 15 ft. -1 M , 
where h = H /(100kms~ 1 Mpc~ 1 ), being H the Hubble constant at current epoch (in 
the paper we adopt h — 1/2). This mass comes from a region of diameter ~ 20/j" 1 Mpc 
and consequently the observations of clusters can provide information on the mass distri- 
bution of the Universe on these scales. Furthermore, since rich clusters are rare objects, 
their properties are expected to be sensitive to the underlying mass density field from 
which they arose. On scales of ~ 20h~ 1 Mpc the density field is still described by the 
linear perturbation theory so the measurement at the present epoch can be related to 
the initial conditions. While the gravitational mass in bound objects is easily computed 
using analytic approaches such as the PS formalism or N-body simulations, it is not so 
easily measured The most robust measurements of the clusters abundance currently rely 
on the clusters temperature function, that is the number density of clusters above a cer- 
tain temperature N(> kT) expressed in units of h 3 Mpc~ 3 . The reasons of this choice 
are several. The integrated temperature does not suffer from projection effects, and com- 
plete samples of clusters may be obtained from all sky X-ray surveys. Clusters are close 
to isothermal, both observationally and in simulations, which makes their temperature 
determination robust and insensitive to numerical resolution or telescope angular resolu- 
tion. The temperature of a cluster depends primarily on the depth of the potential well 
of the dark matter, and the state of equilibrium of the gas. This is in contrast to the 
clusters luminosity, which depends strongly on small scale parameters like clumping and 
core radius. 

The clusters abundance is one of the best observables available for determining the den- 
sity field. In fact, according to the gravitational instability scenario, galaxies and clusters 
form where the density contrast, 5, is large enough so that the surrounding matter can 
separate from the general expansion and collapse. Consequently the abundance of col- 
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lapsed objects depends on the amplitude of the density perturbations. In the CDM model 
these latter follow a gaussian probability distribution and their amplitude on a scale R is 
defined by u(R), the r.m.s. value of 5, which is related to the power spectrum, P(k). In hi- 
erarchical models of structure formation, like CDM, a(R) decreases with increasing scale, 
R, and consequently the density contrast required to form large objects, like clusters of 
galaxies, rarely occurs. The present abundance of clusters is then extremely sensitive to 
small change in the spectrum, P(k). Moreover the rate of clusters evolution is strictly 
connected to the density parameter, flo. Then clusters abundance and its evolution are 
a probe of fl and P(k) and can be used to put some constraints on them. 
Henry and Arnaud (1991) and Edge et al. (1990) have analysed the HEAO-1 A2 all sky 
X-ray survey to obtain temperatures for clusters at a flux limit of F 2 ^iokeV > 3 x 10~ n 
erg/(cm 2 sec). The sample contains the 25 brightest X-ray clusters and is complete for 
more than 90 % at galactic latitude |6| > 20°. By virtue of their brightness the temper- 
atures are easily and well determined. They also constructed temperature distributions 
that has been compared with the CDM predicted X-ray temperature distribution func- 
tion (Bartlett & Silk 1993). Using this data it was shown that galaxy clusters are too 
numerous in the CDM scenario (Bartlett & Silk 1993). 

The other problem of the CDM model that we are addressing is that of the discrepancy 
between the two-point correlation function of clusters of galaxies and the observed one. 
Measurements of galaxy clustering on large scales (Maddox et al 1990a; Efstathiou et 
al. 1990b; Saunders et al. 1991) revealed rms fluctuations of the order of 50% within 
spheres of radius 2Qh~ 1 Mpc. These amplitudes, relative to non-linear clustering on scales 
5/i _1 Mpc, are 2-3 times larger than that predicted by the CDM model. At first the sug- 
gestion of excess power arose from the estimates of the autocorrelation of Abell clusters 
(Hauser & Peebles 1973; Bahcall & Soncira 1983; Klypin & Kopylov 1983). As shown 
by White et al. (1987), the cluster correlations in Standard Cold Dark Matter (SCDM) 
should have an amplitude 2-3 times smaller than these estimates. Although the integrity 
of Abell cluster has been repeatedly called into question (Sutherland 1988; Dekel et al. 
1989; Olivier et al. 1990; Sutherland & Efstathiou 1991) and a new sample of rich clusters 
identified from the APM survey exhibits weaker correlations, marginally consistent with 
the CDM predictions (Dalton et al. 1992), X-ray clusters samples (Lahav et al. 1989) 
have a large correlation length and observations of bright radio galaxies (Peacock 1991; 
Peacock & Nicholson 1991) are also strongly clustered on large scales. 
Finally, as shown in some studies of galaxy clustering on large scales (Maddox et al. 
1990a; Efstathiou et al. 1990b; Saunders et al. 1991), the measured rms fluctuations 
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within spheres of radius 20h~ 1 Mpc have values 2-3 times larger than that predicted by 
the CDM model. In order to solve the quoted discrepancies between CDM model previ- 
sions and observations, alternative models have been introduced. Several authors (Peebles 
1984; Efstathiou et al. 1990a; Turner 1991) have lowered the matter density under the 
critical value (f2 TO < 1) and have added a cosmological constant in order to retain a flat 
Universe (O m + J1a = 1)- The spectrum of the matter density is specified by the transfer 
function, but its shape is affected because the epoch of matter-radiation equality (charac- 
terized by the redshift z eq ) is earlier, l + z eq being increased by a factor l/fi m . Around the 
epoch of redshift z\, where z\ — (Q^/Qm) 1 / 3 — 1, the effect of the cosmological constant 
becomes important and the growth of the density contrast slows down and ceases after 
z\. As a consequence, the normalisation of the transfer function begins to fall, even if its 
shape is retained. In particular Bartlett & Silk (1993) showed that a model with £1 = 0.2 
and tt\ = 0.8 produces a reasonable temperature distribution of clusters, but a higher 
normalization with <r 8 ~ 2, (where as is the rms value of in a sphere of 8h~ 1 Mpc), is 
needed to explain the peculiar velocity field (Efstathiou et al. 1992b). 
Mixed dark matter models (MDM) (Shafi & Stecker 1984; Valdarnini & Bonomctto 1985; 
Schaefer et al. 1989; Holtzman 1989; Schaefer 1991; Schaefer & Shafi 1993; Holtzman & 
Primack 1993) increase the large-scale power because neutrinos free-streaming damps 
the power on small scales. In particular the model with Qhdm — 0.3 simulated by Davis 
et al (1992), gives a good reproduction of the clusters abundance with a normalization 
within the 1 a errors of the COBE. These last models have some difficulty in reproducing 
peculiar velocities (Efstathiou et al. 1992b). 

Alternatively, changing the primeval spectrum, several problems of CDM are solved (Cen 
et al. 1992). Finally it is possible to assume that the threshold for galaxy formation is not 
spatially invariant but weakly modulated (2% — 3% on scales r > 10h~ 1 Mpc) by large 
scale density fluctuations, with the result that the clustering on large-scale is significantly 
increased (Bower et al. 1993). This model follows the spirit of the well known high-peak 
model but differs from it because nonlocal physical processes produce different shapes of 
the mass and galaxy correlation function. 

In any case the solution to the quoted problems till now proposed is related to alternative 
models with more large-scale power than CDM. 

Here, we propose a solution to the problem using the CDM model and taking account of 
the non-radial motions originating from the gravitational interaction of the quadrupole 
moment of the protocluster with the tidal field of the matter of the neighboring proto- 
structures. 
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The plan of this work is the following: in Sect. 2 we introduce the model used to find 
the effects of non-radial motions on the X-ray temperature distribution function and the 
two-point correlation function. In Sect. 3 we compare the results of this model with 
the X-ray temperature distributions given by Henry & Arnaud (1991) and Edge et al. 
(1990) while in Sect. 4 we compare our model with the two-point correlation function 
obtained by Sutherland & Efstathiou (1991), from the analysis of Huchra's et al. (1990) 
deep redshift survey, and with the data points for the Automatic Plate Measuring (APM) 
clusters, computed by Efstathiou et al. (1992a). Sect. 5 is devoted to the conlusions. 



2. The X-ray temperature function 

The PS theory provides an analytical description of the evolution of structures in a 
hierarchical Universe. In this model the linear density field, p(x, t), is an isotropic random 
Gaussian field, the non-linear clumps are identified as over-densities (having a density 
contrast S c ~ 1.68 - Gunn & Gott 1972) in the linear density field, while a mass element 
is incorporated into a non-linear object of mass M when the density field smoothed with 
a top-hat filter of radius Rf, exceeds a threshold S c (M oc Rj)- The comoving number 
density of non-linear objects of mass M to M + dM is given simply by: 

( ' *)<™ - -/*Vr exp (-72) - [ m ) w (1) 

where pb is the mean mass density, a(M) is the rms linear mass overdensity evaluated at 
the epoch when the mass function is desidered and v = g ^ . The redshift dependence 
of Eq. (1) can be obtained remembering that 

_ 5 c (z)D(0) 

a (M)D(z) [ ' 

being D(z) the growth factor of the density perturbation and <r (M) the current value 
of a{M). In Eq. (1) PS introduced arbitrarily a factor of two because dF(M) = 1/2, 
so that only half of the mass in the Universe is accounted for. Bond et al. (1991) showed 
that the "fudge factor" 2 is naturally obtained using the excursion set formalism in the 
sharp fc-space while for general filters (e.g., Gaussian or "top hat") it is not possible 
to obtain an analogous analytical result. As stressed by Yano et al. 1996, the factor of 
2 obtained in the sharp fc-space is correct only if the spatial correlation of the density 
fluctuations is neglected. In spite of the quoted problem, several authors (Efstathiou et 
al. 1988; Brainerd & Villumsen 1992; Laccy & Cole 1994) showed that PS analytic theory 
correctly agrees with N-body simulations. In particular Efstathiou et al. (1988), showed 
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that PS theory correctly agrees with the evolution of the distribution of mass amongst 
groups and clusters of galaxies (multiplicity function). Brainerd & Villumsen (1992) stud- 
ied the CDM halo mass function using a hierarchical particle mesh code. From this last 
work it results that PS formula fits the results of the simulation up to a mass of 10 
times the characteristic ltr fluctuation mass, M*, being M* ~ 10 15 6 _6 ^™ !+3 - ) Mo, where 
b is the bias parameter and ri; is the local slope of the power spectrum. PS theory has 
proven particularly useful in analyzing the number counts and redshift distributions for 
QSOs (Efstathiou & Rees 1988), Lyman a clouds (Bond et al. 1988) and X-ray clusters 
(Cavaliere & Colafrancesco 1988). 

Some difficulties arise when PS theory is compared with observed distributions. To esti- 
mate the multiplicity function of real systems one needs to know the temperature-mass 
(T-M) relation in order to trasform the mass distribution into the temperature distri- 
bution. Theoretical uncertainty arises in this transformation because the exact relation 
between the mass appearing in the PS expression and the temperature of the intracluster 
gas is unknown. Under the standard assumption of the Intra-Cluster (IC) gas in hydro- 
static equilibrium with the potential well of a spherically simmetric, virialized cluster, 
the IC gas temperature-mass relation is easily obtained by applying the virial theorem 
and for a flat matter-dominated Universe we have that (Evrard 1990, Evrard et al. 1996, 
Evrard 1997, Bartlett 1997): 
M 



T=(6.4*V3 fce ^__j (1 + z) (3 ) 

The assumptions of perfect hydrostatic equilibrium and virialization arc in reality not 
completely satisfied in the case of clusters. Clusters profile may depart from isothermality, 
with sligth temperature gradients throughout the cluster. The X-ray weighted temper- 
ature can be slightly different from the mean mass weighted virial temperature. In any 
case the scatter in the T-M relation given by Eq. (3) is of the order of ~ 10% (Evrard 
1991). The mass variance present in Eq. (1) can be obtained once a spectrum, P(k), is 
fixed: 

a 2 {M) = —z \ dkk 2 P{k)W 2 {kR) (4) 
27r Jo 

where W(kR) is a top-hat smoothing function: 

W(kR) = -^-^ (sinfci?- kRcoskR) (5) 
and the power spectrum P(k) = Ak n T 2 (k) is fixed giving the transfer function T(k) : 
T(k) = [In (1 + 4.164fc)] 2 • (192.9 + 1340fc + 

+ 1.599 • 10 5 fc 2 + 1.78 • 10 5 fc 3 + 3.995 • 10 6 fc 4 )~ 1/2 (6) 
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(Ryden & Gunn 1987; BBKS) and A is the normalizing constant. The normalization 
of the spectrum may be obtained in several ways. One possibility is to normalize it to 
COBE scales using the cosmic microwave anisotropy quadrupolc Q rms -ps — l7(iK. This 
corresponds to a 8 = 0.95 ± 0.2 (Smoot et al. 1992; Liddlc & Lyth 1993). More recent 
determinations give a% = 1-33 if we use the BBKS spectrum, while a$ — 1.22 if we use the 
spectrum by Bond & Efstathiou (1984) (Klypin et al. 1997). Another way of fixing the 
normalization is via the abundance of clusters giving a$ = 0.5 -j- 0.6 (Pen, 1997; Oukbir 
1996; Bartlett 1997). Normalisation on scales from 10 to 50Mpc obtained from QDOT 
(Kaiser et al. 1991) and POTENT (Dekel et al. 1992) requires that a 8 is in the range 
0.7-j- 1.1, which is compatible with COBE normalisation while the observations of the 
pairwise velocity dispersion of galaxies on scales r < 3Mpc seem to require a 8 < 0.5. Our 
normalization, erg = 1, is intermediate between that suggested by clusters abundance and 
that of COBE. 

As shown by Bartlett & Silk (1993) the X-ray distribution function obtained using a stan- 
dard CDM spectrum over-produces the clusters abundances data obtained by Henry & 
Arnaud (1991) and Edge et al. (1990). The discrepancy can be reduced taking into account 
the non-radial motions that originate when a cluster reaches the non-linear regime. In fact, 
the PS temperature distribution requires the specification of 5 C and the temperature-mass 
relation T-M. The presence of non-radial motions changes both S c and the T-M relation. 
Barrow & Silk (1981), Szalay & Silk 1983 and Peebles 1990 assumed that non-radial 
motions would be expected (within a developing protocluster) due to the tidal interac- 
tion of the irregular mass distribution around the protocluster with the neighbouring 
protoclusters. The kinetic energy of these non-radial motions inhibits the collapse of the 
protocluster enabling it to reach statistical equilibrium before the final collapse (Davis 
& Peebles 1977; Peebles 1990). The role of non-radial motions has been also pointed 
out by Antonuccio & Colafrancesco (1995). After deriving the conditional probability 
distribution f p k(v\v) of the peculiar velocity around a peak of a Gaussian density fiery 
non-radial motions. In these regions the fate of the infalling material could be influenced 
by the amount of tangential velocity relative to the radial one. 

This can be shown writing the equation of motion of a spherically symmetric mass dis- 
tribution with density n(r): 

| nK) + JU^) + ( 2{v 2 } _ {v 2 }) 2 + n(r)|K) = (7) 

where (v r ) and (v#) are, respectively, the mean radial and tangential streaming velocity. 
Eq. (7) shows that high tangential velocity dispersion ((v#) > 2{v%)) may alter the infall 
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pattern. The expected delay in the collapse of a perturbation may be calculated using a 
model due to Peebles (Peebles 1993). 

Let's consider an ensemble of gravitationally growing mass concentrations and suppose 
that the material in each system collects within the same potential well with inward 
pointing acceleration given by g(r). We indicate with dP = /(£, rv r , t)dLdv r dr the prob- 
ability that a particle can be found in the proper radius range r, r + dr, in the radial 
velocity range v r = r, v r + dv r and with angular momentum L — rvg in the range dL. 
The radial acceleration of the particle is: 



dv r L 2 (r, v) 
~dl = M 2 r 3 



(r) 



(8) 



Eq. (8) can be derived from a potential and then from Liouville's theorem it follows that 
the distribution function, /, satisfies the collisionless Boltzmann equation: 



d_l +v d_l + dl 

dt dr dv r 







(9) 



Using Gunn & Gott's (1972) notation we write the proper radius of a shell in terms of 
the expansion parameter, a(r i} t), where is the initial radius: 



r(n,t) = na(ri,t) 



(10) 



Remembering that M = ^-pfa, t)a 3 (ri, t)rf, that |^ = p C i, with p ci and Hi respectively 
the critical mass density and the Hubble constant at the time tj, and assuming that no 
shell crossing occurs so that the total mass inside each shell remains constant, (p(rj, t) = 

d 6 (n,t) 



p,(r "*' ) ) Eq. (8) may be written as: 



fa = Hf(l + S) 

dt 2 ~ ~ " 



2a 2 



+ 



AG 2 L 2 



Hf(l + S) 2 rj°a 3 
Integrating Eq. (11) we have: 



(11) 



where 6 = Pi p " 

Pci 



8G 2 L 2 



da-2C 



where C is the binding energy of the shell. Integrating once more we have: 

da 



I 

Jo 



H 2 



1+5 



l+S 



+ r 



8G 2 L 2 



(12) 



(13) 



Using Eqs (12) and (13) it is possible to find the linear over-density at the turn-around 
epoch, tta- In fact solving Eq. (13), for some epoch of interest, we may obtain the ex- 
pansion parameter of the turn-around epoch. This is related to the binding energy of 
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the shell containing mass M by Eq. (12) with ^| = 0. In turn the binding energy ol a 
growing mode solution is uniquely given by the linear overdensity, Si, at time ij. From 
this overdensity, using the linear theory, we may obtain that ol the turn-around epoch 
and then that of the collapse which is given by: 

8G 2 f a ^L 2 -da] 



5 c {v) 



1 



(14) 



W o H$r} S(l + <5) 2 Jo a? 
where S co = 1.68 is the critical threshold for a spherical model, while H and il are 
respectively the Hubble constant and the density parameter at the current epoch to- 
Filtering the spectrum on clusters scales, Rf — 3h~ 1 Mpc, we obtained the total specific 
angular momentum, h{r,v) = L{r,v)jM sh , acquired during expansion, integrating the 
torque over time (Ryden 1988 - Eq. 36): 



r o t s o r (i-cosg) J Mfrdd 

[ ^8M sh Jo (<? - sin ) 4/3 h (#) - h W *f 



where t S a and S are respectively the torque, the mean overdensity and the mean 
overdensity within a sphere of radius r at the current epoch to. The functions fi("&), 
/ 2 (i?) are given by Ryden (1988 - Eq. 31): 

f^O) = 16-16cos6» + sin 2 6»-96»sin0 (16) 

f 2 (9) = 12 - 12 cos 6 + 3 sin 2 9 - 99 sin 6 (17) 

where 9 is a a parameter connected to the time, t, through the following equation: 

t=^t S~ 3/2 (9 -sin 9) (18) 

The mean overdensity within a sphere of radius r , S(r), is given by: 

S(r,v) = 4 / dxx 2 S(x) (19) 
r Jo 

The mass dependence of the threshold parameter, 5 c {v), can be found as follows: we 
calculate the binding radius, r^, of the shell using Hoffmann & Shaham's criterion (1985): 



T c (r, v) < to (20) 

where T c (r, v) is the calculated time of collapse of a shell and t Q is the Hubble time. We 
find a relation between v and M through the equation M = 4ir Pbi*l / '3 . We so obtain 
8 c {y(M)). In Fig. 1 we show the variation of the threshold parameter, S C (M), with the 
mass M. Non-radial motions influence the value of S c which increases for peaks of low mass 
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Fig. 1. The threshold 5 C as a function of the mass M, for a CDM spectrum (fio = 1, h = 1/2) 
with Rf — 3/i _1 Mpc, taking account of non-radial motions. 

and remaines unchanged for high mass peaks. As a consequence, the structure formation 
by low mass peaks is inhibited. In other words, in agreement with the cooperative galaxy 
formation theory (Bower et al. 1993), structures form more easily in over-populated 
regions. As we previously told, the cooperative galaxy formation is able to reconcile the 
CDM model with the APM correlations by assuming the threshold for galaxy formation 
to be modulated by large-scale density fluctuations rather than to be spatially invariant. 
But there exists some difficulty in finding a physical mechanism able to produce the 
modulation. In our model this mechanism is linked to non-radial motions. 
To get the temperature distribution it is necessary to know the temperature-mass relation. 
This can be obtained using the virial theorem, energy conservation and using Eq. (12). 
From the virial theorem we may write: 

GM [ r ' ss L 2 
while from the energy conservation: 
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-(K) 



GM 

r eff 



+ 



• r *ff 



L 1 



M 2 r 3 



dr 



gm r 

r ta Jo 



M 2 r 3 



dr 



Eq. (21) and Eq.(22) can be solved for r e // and < K >. We finally have that: 



T = (6AkeV) 



(- 



M-h 



V1O 15 M 



2/3 



r L'dr 
H 1 ^ 



where r\ is a parameter given by r\ — r ta /xi, being r ta the radius of the turn-around 
epoch, while x\ is defined by the relation M = 4irpi,Xi/3 and ip — r e ff/r ta where r e // is 
the time-averaged radius of a mass shell. Eq. (23) was normalised to agree with Evrard's 
(1990) simulations for L = 0. 



(22) 



(23) 



3. Non radial motions and the X-ray temperature function. 

The new T-M relation is Eq. (23) which differs from Eq. (3) for the presence of the term: 
^ Jo ^24) 

( G 2;3£o M 2)l/3 

This last term changes the dependence of the temperature on the mass, M, in the T-M 
relation. Moreover the new T-M relation depends on the angular momentum, L, origi- 
nating from the gravitational interaction of the quadrupole moment of the protocluster 
with the tidal field of the matter of the neighboring protostructures. In Fig. 2 the X-ray 
temperature distribution, derived using a CDM model with Clo = 1, h = 1/2 and taking 
into account non-radial motions, is compared with Henry & Arnaud (1991) and Edge et 
al. (1990) data and with a pure CDM model with fl = l,h= 1/2. As shown the CDM 
model that does not take account of the non-radial motions over-produces the clusters 
abundance. The introduction of non-radial motions gives a more careful description of 
the experimental data. As we have seen the X-ray temperature distribution function ob- 
tained taking account of non-radial motions is different from that of a pure CDM model 
for two reasons: 

1) the variation of the threshold, S c , with mass, M. This is due to the change in the en- 
ergetics of the collapse model produced by the introduction of another potential energy 
term { L ^2 V \ ) i n t ne equation decsribing the collapse (see Eq.(8)): 

2) the modification of the T-M relation produced by the alteration of the partition of 
energy in virial equilibrium. 

For values of mass M = 0.5M Q the difference between the two theoretical lines in Fig. 
2 is due to the first factor for ~ 59% and this value increases with increasing mass. The 
uncertainty in our model fundamentally comes from the uncertainty of the T-M relation 
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Fig. 2. X-ray temperature distribution function. The dashed line gives the temperature function 
for a pure CDM model (fio = 1, h = 1/2), with Rf = 3h~ 1 Mpc. The solid line is the same 
distribution but now taking account of non-radial motions. The data are obtained by Edge et 
al. 1990 (dots), and Henry & Arnaud 1991 (open squares) 

whose value has been previously quoted. 

Somebody may object that the effect here described has not been seen in some hy- 
drodynamic simulation (Evrard & Crone 1992). The answer to this objection can be 
given remembering a similar problem of the previrialization conjecture (Davis & Peebles 
1977; Peebles 1990), (supposing that initial asphericities and tidal interactions between 
neighboring density fluctuations induce significant non-radial motions, which oppose the 
collapse) on which our model is fundamentally based. It is known that while some N-body 
simulations (Villumscn & Davis 1986; Peebles 1990) appear to reproduce this effect, other 
simulations (for example Evrard & Crone 1992) do not. An answer to this controversy 
was given by Lokas et al. 1996. The problem is connected to the spectral index n used 
in the simulations. The "previrialization" is seen only for n > —1. While Peebles (1990) 
used simulations with n = 0, Evrard & Crone (1992) assumed n = —1. Excluding this 
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particular case, generally the whole properties of clusters, such as their optical and X-ray 
luminosity functions, or their velocity and temperature distribution functions, are diffi- 
cult to address directly in numerical simulations because the size of the box must be very 
large in order to contain a sufficient number of clusters. Then the analytical approach 
remains an effective alternative. 



4. Non-radial motions and the clusters correlation function 

As previously seen in Sect. 2, in the PS theory the comoving number density of non-linear 
objects of mass M to M + dM is simply obtained by differentiating with respect to mass 
the integral from 5 C to infinity of the probability distribution for fluctuations given by: 

p[8(M)\ = ^ exph5(M) 2 /2 ( r(M) 2 ] (25) 

In an exactly analogous way, the probability p(M\, M 2 , r) per unit volume per unit masses 
of finding two collapsed objects of mass Mi and M 2 separated by a distance r is obtained 
by integrating both variables of the bivariant Gaussian distribution in 5 (Mi) and 6(M 2 ), 
with correlation £ p (r) from S c to infinity and then taking the partial derivatives with 
respect to both masses. The correlation function for collapsed objects is simply obtained 
from: 

£,mm = p(Mi,M 2 , r)/p(Mi)p(M 2 ) - 1 (26) 

which for equal masses and on scales bigger than that of the turn-around (£ p <C 1) is 
(Kashlinsky 1987): 

Hmm = [5 2 Ja(Mf%(r) (27) 

where £ p is the correlation function of the matter density distribution when the den- 
sity fluctuations had small amplitude. Eq. (27) shows that the correlation of collapsed 
objects may be enhanced with respect to the underlying mass fluctuations. This con- 
dition is usually described by the bias parameter b which is sometimes defined as 
lZMM(r)/Ur)} 1/2 =5c/a(M). 

Studies of clustering on scales > 10h~ 1 Mpc have shown that the correlation function 
given in Eq. (27) is different from that obtained from observations. The most compelling 
data are angular correlation functions for the APM survey. These decline much less 
rapidly on large scales than the CDM prediction (Maddox et al. 1990a). This discrepancy 
can be reduced taking into account the non-radial motions that originate when a cluster 
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reaches the non-linear regime. In fact, the calculation of the correlation function requires 
the specification of 5 C which is modified by non-radial motions. 
The result of our calculation is showed in Fig. 3. 

The observational data to which we compare the calculated correlation function are those 
obtained by Sutherland & Efstathiou (1991) from the analysis of Huchra's et al. (1990) 
deep redshift survey as discussed in Geller & Huchra (1988) and the data points for the 
APM clusters computed by Efstathiou et al. (1992a). The deep-cluster redshift survey of 




r (Mpc) 



Fig. 3. Clusters of galaxies correlation function. The solid line gives the correlation function for 
a pure CDM model, with Rf — 3hr 1 Mpc. The dashed line is the same distribution but now takes 
account of non-radial motions. The observational data refer to the two-point correlation function 
obtained by Sutherland & Efstathiou (1991) [filled exagons) from the analysis of Huchra's et 
al. (1990) deep redshift survey and with the data points for the APM clusters computed by 
Efstathiou et al. (1992a) [dashed errorbars). 

Huchra et al. (1990) consists of the 145 Abell clusters with R > 0, D < 6 in the area 
W h < a < 15'\ 58° < S < 78°. The APM sample consists of 240 clusters with APM 
richness R > 20 and photometrically derived redshifts z x < 0.1. The survey covers a solid 
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angle w = l.Ssr in the region of sky 21 h < a < 5 h , -72.5° < S < -17.5° of the APM 
galaxy survey. Both in Sutherland & Efstathiou (1991) and Efstathiou et al. (1992a) the 
correlation function was calculated from the samples using the estimator: 

W= F W L - 1 (28) 

where N cc and N cr are the numbers of cluster pairs and cluster-random pairs having 
redshift-space separations in the range (r — dr/2,r + dr/2). The random points are gen- 
erated within the sample volume with a mean density F times that of the clusters and 
with a redshift distribution derived from a smoothed distribution of the redshifts for the 
cluster sample. The correlation function obtained was finally corrected for line-of-sight 
anisotropies. This correction is necessary because the clustering of Abell clusters is highly 
anisotropic in redshift space, providing evidence that the Abell catalogues, which were 
built by scanning photographic plates by eye, are affected by incompleteness on the plane 
of the sky which enhances the clustering amplitude measured in three dimensions (Suther- 
land & Efstathiou 1991; Efstathiou et al. 1992a). In particular Efstathiou et al. (1992a) 
presented a clear example of incompleteness in the Abell R > catalogue comparing the 
machine-based APM survey with the Postman, Hucra & Geller 1992 redshift survey of 
Abell clusters showing that even if the clusters in these surveys have comparable space 
densities the red-shift space correlation function of the APM sample is isotropic on large 
scales while the correlation function for the Abell clusters is highly anisotropic. In any 
case, after this correction the correlation function of Abell clusters agrees extremely well 
with that of APM clusters (Efstathiou et al. 1992a). We use two samples because: 
a) we wanted to compare our model for the two-point correlation function with observa- 
tional data and at the same time with the result obtained by Borgani (1990) who studied 
the effect of particular thresholds (erfc-threshold and Gaussian-threshold) on the corre- 
lation properties of clusters of galaxies and compared his result with that by Sutherland 
& Efstathiou (1991). For this reason we used this last sample; 

2) we used the APM galaxy survey because the uniformity of the APM magnitudes, the 
low obscuration in the APM survey area and the use, for its construction, of a computer 
cluster-finding algorithm (see Maddox et al. 1990a, b) further reduce the possibility of 
spurious clustering on the plane of the sky. 

The comparison between the correlation function and the quoted data, displayed in Fig. 
3, shows that there is an evident discrepancy between pure CDM previsions and exper- 
imental data. The CDM model seems to have trouble in re-producing the behaviour of 
the data. In fact, the predicted two-point cluster function is too steep and rapidly goes 



17 



nearly to zero for r ~ 30/i _1 Mpc, while the data show no significant anticorrelation up to 
r ~ 60h~ 1 Mpc (see Borgani 1990). This is a direct effect due to the non-scale invariance 
of the CDM spectrum (Primack & Blumcnthal 1983). The introduction of non-radial mo- 
tions gives a more accurate description of the experimental data showing that physical 
effects cannot be ignored in the study of the formation of cosmic structures. The result 
obtained is in agreement with that obtained by Borgani (1990) who showed how the 
introduction of smooth thresholds (similar to that obtained in our model) leads to two- 
point correlation functions in a systematically agreement with the data. Before going on 
we want to remember that the threshold functions are strictly connected to the concept 
of bias. In fact, according to the biased theory of galaxy formation, observable objects of 
mass ~ M arise from fluctuations of the density field, filtered on a scale R, rising over 
a global threshold, 8 > 8 C — v t o, where a is the rms value of 8 and v t is the threshold 
height. The number density of objects, n p k, that forms from peaks of density of height v 
can be written following BBKS (1986) in the form: 

rOO 

n pk = / t(-)N pk (v)dv (29) 
Jo v t 

where t( j^) is the threshold function, v t the threshold height and N p kdv the differential 
number density of peaks (see BBKS 1986 Eq. 4.3). The threshold level v t is defined so 
that the probability of a peak becoming an observable object is 1/2 when v — v t . In the 
sharp threshold case the selection function, is a Hcaviside function t(-jf-) = 6(y — Vt). As 
previously quoted the threshold function is connected to the bias coefficient of a class of 
objects by (BBKS): 

6(^ = 1^ + 1 (30) 

(To 

where < v > is: 

70 



< 



Jo 



t{-)N pk {v)dv (31) 



1 — 7 : 

while, 7 and •& are given in BBKS (1986) (respectively Eq. 4.6 a; Eq. 6.14). 
While in a 9 threshold scheme, fluctuations below 8 C have zero probability to develop 
an observable object and fluctuations above 8 C have zero probability not to develop an 
object, the situation is totally different when an erfc-threshold, as that introduced by 
Borgani (1990), is used. In this case objects can also be formed from fluctuations below 
8 C and there is a non-zero probability for fluctuations above 8 C to be sterile. 
According with Borgani (1990), the erfc-threshold can be related to non-sphericity effects 
during the gravitational growing process. In fact, while in the spherical limit it is possible 



18 



to find a precise relation between the time, t elapsed from the turning around and the 
correspondent density contrast, 5,{t oc S~ 3 / 2 ) (Gunn & Gott 1972), when non-sphericity 
is introduced it is no longer possible to univocally relate the primeval density contrast 
and the evolutionary stage of a fluctuation because while in the spherical model an object 
is characterized by 5 > S c the non-sphericity produces a spread around the typical value 
5 C . The erfc-threshold is a way to correct the quoted t — 5 relation. The erfc-threshold is 
also linked to non-radial motions. In fact as we previously told, the tidal interaction of the 
irregular mass distribution within and around the protocluster, present in hierarchical 
models, gives rise to non-radial motions. We also told that the erfc-threshold is related to 
non-sphericity effects that in turn are responsible for the origin of the quoted non-radial 
motions. 

Borgani's (1990) model, like our, is characterized by a non-6* threshold. As we showed in a 
previous paper, Del Popolo & Gambera 1998a, one of the effects of non-radial motions is 
that the threshold function differs from a Heaviside function (sharp threshold), (see Fig. 
7 of Del Popolo & Gambera 1998a). In this last paper the threshold function is defined 
as: 



gives the probability that the peak overdensity is different from the average, in a Gaussian 
density field. As displayed, the integrand is evaluated at a radius rut which is the typical 
radius of the object that we are selecting. Moreover, the threshold function t(v) depends 
on the critical overdensity threshold for the collapse, <5 C , which is not constant as in a 
spherical model (due to the presence, in our analysis, of non-radial motions that delay 
the collapse of the proto-cluster) (see Eq. (14)). The fundamental difference between 
our and Borgani's approach is that our threshold function is physically motivated: it is 
simply obtained from the assumptions of a Gaussian density field and taking account 
of non-radial motions. Borgani's threshold functions (erfc and Gaussian threshold) are 
ad-hoc introduced in order to reduce the discrepancy between the observed and the CDM 
predicted two-point correlation functions of clusters of galaxies. The connection with the 
quoted non-sphericity effects, even if logical and in agreement with our results, is only a 
posteriori tentative to justify the choice made. 




(32) 



where the function 




(33) 
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5. Conclusions 

In these last years many authors (Bahcall & Soneira 1983; White et al. 1987; Maddox 
et al. 1990a; Saunders et al. 1991; Peacock 1991; Bartlett & Silk 1993) have shown the 
existence of a strong discrepancy between the observed X-ray temperature distribution 
function of clusters and that predicted by a CDM model and the observed two-point 
correlation function of clusters and that predicted by the CDM model. To reduce these 
discrepancies several alternative models have been introduced but no model has consid- 
ered the role of the non-radial motions. Here we have shown how non-radial motions may 
reduce both the two quoted discrepancies. To this aim we calculated the variation in the 
treshold parameter, 5 C , as a function of the mass M, and that of the temperature-mass 
relation, produced by the presence of non-radial motions in the outskirts of clusters of 
galaxies. We compared Edge et al. (1990) and Henry & Arnaud (1991) data with the dis- 
tribution function of X-ray temperature, calculated using PS (1974) theory and Evrard's 
(1990) prescriptions for the mass-temperature relation. We found that the model pro- 
duces a reasonable clusters temperature distribution (see Fig. 2). We also used 5 C (M) to 
calculate the two-point correlation of clusters of galaxies and we compared it both with 
that obtained by Sutherland & Efstathiou (1991) from the analysis of Huchra's et al. 
(1990) deep redshift survey as discussed in Gcller & Huchra (1988) and with the data 
points for the APM clusters computed by Efstathiou et al. (1992a). Our results (see Fig. 
3) show how non-radial motions change the correlation length of the correlation function 
making it less steep than that obtained from a pure CDM model where the non-radial 
motions are not considered. 
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